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Abstract The quantum dynamics of a spin-1/2 charged par¬ 
ticle in the presence of magnetic field is analyzed for the 
general case where scalar and vector couplings are consid¬ 
ered. The energy spectra are explicitly computed for differ¬ 
ent physical situations, as well as their dependencies on the 
magnetic field strength, spin projection parameter and vec¬ 
tor and scalar coupling constants. 

1 Introduction 

The study of relativistic quantum systems under the influ¬ 
ence of magnetic field and scalar potentials has attracted 
attention of researchers in various branches of physics. It 
is well known that these potentials can be inserted into the 
Dirac equation 

[)37-p + j3M] v/(r) =£'i/(r), (1) 

through three usual substitutions, known as minimal, vec¬ 
tor and scalar couplings, whose representation is denoted, 
respectively, by 

p^p-eA, (2) 

E^E-V{r), (3) 

M^M + S{r). (4) 

a variety of relativistic and nonrelativistic effects can be stud¬ 
ied. Those couplings above differ in the manner how they are 
inserted into the Dirac equation [1]. The minimal coupling 
(2) is useful for studying the dynamics of a spin-1/2 charged 
particle in a magnetic field. For example, using this model, 
we can study Landau levels [2], Aharonov-Bohm effect [3], 
Hall effect [4], and other effects associated with magnetic 
field. 
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It is known that the prescription (3) acts differently on 
electron and positron states, respectively, and the eigenvalue 
spectrum of the particle is not symmetric. In this case, bound 
states exist for only one of the two kinds of particles. In other 
words, we can say that, for vector coupling, the potential 
couples to the charge. In the context of the Dirac equation, 
this coupling has been used, for example, to study the influ¬ 
ence of a harmonic oscillator on the Aharonov-Casher prob¬ 
lem [5], the Aharonov Bohm effect for a spin-1/2 particle in 
the case that a 1/r potential is present [6], effects of non¬ 
gauge potentials on the spin-1/2 Aharonoh-Bohm problem 
[7], quasiclassical theory of the Dirac equation with applica¬ 
tions in the physics of heavy-light mesons [8] and confining 
potentials with pure vector coupling [9]. In the Schrodinger 
theory. It also has important applications, such as, in the dy¬ 
namics of an electron in a two-dimensional quantum ring 
[10, 11], quantum particles constrained to move on a coni¬ 
cal surface [12] and effect of singular potentials on the har¬ 
monic oscillator [13]. 

In the case of the scalar coupling (4), it is added to the 
mass term of the Dirac equation and, therefore, it can be in¬ 
terpreted as an effective, position-dependent mass and, fur¬ 
thermore, it also acts equally on particles and antiparticles. 
This coupling has been used, for example, to obtain an ex¬ 
act solution of the Dirac equation for a charged particle with 
position-dependent mass in the Coulomb held [14], to study 
the relativistic quantum dynamics of a charged particle in 
cosmic string spacetime [15, 16], scattering of a fermion 
in the background of a smooth step potential with a gen¬ 
eral mixing of vector and scalar Lorentz structures with the 
scalar coupling stronger than or equal to the vector cou¬ 
pling [17], inclusion of the generalized Hulthen potential 
in the case of the smooth step mass distribution [18] and 
extension of PT-symmetric quantum mechanics [19]. The 
coupling (4) also has important applications in nonrelativis¬ 
tic quantum mechanics. The Schrodinger equation with a 
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position-dependent mass has attracted a lot of attentions due 
to a wide range of applications in various areas of mate¬ 
rial science and condensed matter physics. For example, to 
study the dynamics of an one-dimensional harmonic oscil¬ 
lator [20], derivation of the Shannon entropy for a particle 
with a nonuniform soli tonic mass density [ 21 ], context of 
displacement operator for quantum systems [22, 23], use of 
instantaneous Galilean invariance to derive the expression 
for the Hamiltonian of an electron [24], determination of 
some potential functions for exactly solvable nonrelativis- 
tic problems [25] and Hermitian, rotationally invariant one- 
band Schrodinger Hamiltonian [26]. 

The case in which the couplings are composed by a vec¬ 
tor (3) and a scalar (4) potentials, with S = V (S = —V), are 
usually pointed out as necessary condition for occurrence of 
exact spin (pseudospin) symmetry. It is known that the spin 
and pseudospin symmetries are SU(2) symmetries of a Dirac 
Hamiltonian with vector and scalar potentials. The pseu¬ 
dospin symmetry was introduced in nuclear physics many 
years ago [27, 28] to account for the degeneracies of or¬ 
bital in single-particle spectra. Also, it is known that the spin 
symmetry occurs in the spectrum of a meson with one heavy 
quark [29] and anti-nucleon bound in a nucleus [30], and the 
pseudospin symmetry occurs in the spectrum of nuclei [31]. 

In this work, we study the quantum dynamics of a spin- 
1/2 charged particle in the presence of magnetic field with 
scalar and vector couplings. This system has been consid¬ 
ered in Ref. [32]. The difference between our approach and 
the previous one is that, here, we solve the problem in a rig¬ 
orous way taking into account other questions which have 
not been analyzed by the authors. For example, absence of 
the term which depends explicitly on the spin in the equa¬ 
tion of motion. Since we are considering the dynamics of a 
particle with spin, such a term can not be neglected in equa¬ 
tion of motion [33]. Moreover, as the authors make a con¬ 
nection with the Aharonov-Bohm problem, the presence of 
this term has important implications on the physical quanti¬ 
ties of interest, such as energy eigenvalues, scattering matrix 
and phase shift (see Ref. [34] for more details). By taking 
into account the term that depends explicitly on the spin into 
Pauli equation, we address the system in connection with the 
spin-1/2 Aharonov-Bohm problem [35] and analyze the fol¬ 
lowing questions: (a) the existence of isolated solutions to 
the first order equation Dirac and (b) general dynamics in 
all space, including the r = 0 region. We use the self-adjoint 
extension method to determine the most relevant physical 
quantities, such as energy spectrum and wave functions by 
applying boundary conditions allowed by the system. 

The paper is organized as follows. In section 2, we con¬ 
sider the Dirac equation in (2 -f 1) dimensions with minimal, 
scalar and vector couplings, and derive the set of first or¬ 
der differential equations. These equations are useful to in¬ 
vestigate possible isolated solutions to the problem. In sec¬ 


tion 3, we solve the first order Dirac equation in connec¬ 
tion with the Aharanov-Bohm problem and scalar and vec¬ 
tor couplings. We found that, for certain values assumed by 
the physical parameters of the system, isolated solutions ex¬ 
ist and discuss the limits of validity of them. In section 4, 
we derive the Pauli equation and study the dynamics of the 
system taking into account exact symmetry spin and pseu¬ 
dospin limits. In section 5, we briefly discuss some concepts 
of the self-adjoint extension method and specify the bound¬ 
ary conditions at the origin which will be used. In section 6 , 
expressions for the energy eigenvalues and wave functions 
are determined for both symmetry limits and compare them 
with the results of Ref. [32]. We verify that the presence of 
the spin element in the equation of motion introduces a cor¬ 
rection in the expressions for the bound state energy eigen¬ 
values. In section 7, we present our concluding remarks. 


2 Equation of motion 


We begin with the Dirac equation (1) in (2-f 1) dimensions 
in polar coordinates {Ti = c = 1) 


{l5r-n + l5[M + S (r)]} V/ (r) = [£ - y (r)] y/ (r), (5) 


where n = (nr, Kif,) = {—idr, —id(p jr — eA^), r = (r, (p) and 
y/ is a two-component spinor. The 7 matrices in Eq. (5) are 
given in terms of the Pauli matrices as [36] 


PY — Gi cos (p -f 5(72 sin <p 


( 0 \ 




= —(7l COS ^ -f S(72 sin ^ 



/ 0 \ 


(6) 

(7) 

( 8 ) 


where s is twice the spin value, with s = 4-1 for spin “up” 
and s = —1 for spin “down”. Equation (5) can be written 
more explicitly as 

[n,- -isn^]\if2 = [E-M-Z (r)] y/i , (9) 

e+“<P y/i = [E+M-A (r)] v/ 2 , (10) 

where E{r) =V (r) + S (r) and A (r) =V (r) — S(r). 

If one adopts the following decomposition 





( 11 ) 


with m+l/2 = ±1/2,±3/2,..., with m G Z, and inserting 
this into Eqs. (9) and (10), we obtain 


d i(m±i) 
- 

dr r 



-esA^ 

-esAff, 


g{r) = [E-M-E{r)]f{r),(U) 

f{r) = [E+M-A{r)]g{r).{n) 
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Note that the above equations are coupled. However, if 
£ (r) or A (r) is made zero in any of the equations, we can 
uncouple them easily. We will see below that this result in 
important physical consequences to the physical system in 
question. 


3 Isolated solutions for the Dirac equation of motion 


In this section, we investigate the existence of isolated solu¬ 
tions in the quantum motion of a fermionic massive charged 
particle in (2-1-1) dimensions. This is accomplished by con¬ 
sidering the particle at rest, i.e., E = ±M, directly in the 
first order equations in Eqs. (12) and (13). Such a solution 
is known to be excluded from the Sturm-Liouville problem, 
and have been investigated under diverse perspectives in the 
latest years [37-44]. We are seeking for bound-state solu¬ 
tions subjected to the normalization condition 


^ {\fm{r)\^ + \gm(.r)\^)rdr=l. 


(14) 


In order to determine the isolated bound-state solutions, 
we consider Z (r) = 0 in Eq. (12), so that, for £ = M, we 
can write 

d s{m + s) 


dr 


— seA 


gm(r)=0, 


d sm 

- — -I- seA 

dr r 


fm{r)=2{M-V {r))gm{r), 


whose general solutions are 

g„{r) = 


Ur) = 

where a+ and b+ are constants, and I{r) is given by 
I{r) = J dr [2M - 2V (r)] 


(15) 

(16) 


(17) 

(18) 


(19) 


which, for a given V (r) and A^ (r), it can be expressed in 
terms of the upper incomplete Gamma function [45] 


"{a,x) = J *e ’dt, 91(a) >0. 


( 20 ) 


Let us now analyze solutions for E = —M, and consider 
A{r) =0 in Eq. (13). Eor this case, we write 


d 

dr 


; (in -f s) 


-esA^ (r) 


g(r) = -2[M + V(r)]f(r), 


d sm 

- — H- esAa (r) 

dr r 


fir)=0, 


whose general solution is 

Ur) = 

gm(r) = 


( 21 ) 

( 22 ) 


(23) 

(24) 


where 

H(r) = J dr [2M - 2V (r)] (25) 

Now, lets us consider the particular case where the parti¬ 
cle moves in a constant magnetic field and in the presence of 
Aharonov-Bohm effect. The vector potential in the Coulomb 
gauge is 


A = Ai -f Az, 

(26) 

with 


Bor ^ (j) ^ 

Ai = —(p, A2 = -(p, 

2 r 

(27) 


where Bq is the magnetic field magnitude and <j) is the flux 
parameter. The potentials in Eq. (26) both provide one mag¬ 
netic field perpendicular to the plane (r, <p), namely 


B = Bi-fB2, (28) 


with 


Bi=VxAi=Boz, (29) 

Bz = VxA2 = ())^z, (30) 

r 

where Biis an external magnetic field and Bz is the magnetic 
field due to a solenoid. If the solenoid is extremely long, the 
field inside is uniform, and the field outside is zero. How¬ 
ever, in a general dynamics, the particle is allowed to access 
the r = 0 region. In this region, the magnetic field is non¬ 
null. If the radius of the solenoid is ro « 0, then the relevant 
magnetic field is Bz ~ 5(r) as in Eq. (30). This situation 
has not been accomplished in Ref. [32], which is crucial 
to give meaning to the term that explicitly depends of the 
spin, namely, the Pauli term that appearing in the second or¬ 
der differential equation. This issue will be considered later 
when we treat solutions for the case E ^ -EM. 

Using Eqs. (29) and (30), we have 


J A(pdr = 


bU 

4 


E <j) In r. 


(31) 


If 0 > 0 and Bq > 0, for E = M, we have bound-state solu¬ 
tions only in the following cases: 


fm(r) 

gm{r) 

fm(r) 

gm{r) 


0 + ’ U+=0, ’ 


(32) 


^ la fjn-X-l -Sr^ ) ■* 

ij“+" ^ ’U+=M = y = o, 


and for E = —M, 
'0 


fm{r) 

gm(r) 

fm(r) 

gm(r) 


I J \ a- =0, 


0J ^ ’ ]b+ =M = V =0, 


(33) 

(34) 

(35) 
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where 5 = eBojA and X = e^. Note that the above results 
are independent of the values of s, m and X to ensure a 
bound state. This is because the function predomi¬ 
nates over the polynomials and If we consider 

Bo —V {f) = 0, which leads to the usual Aharonov-Bohm 
effect, the solutions for £ =M reads 

(36) 

f,n{r) = [b+-a+I{r)]r<’«-^\ (37) 

where 


4.1 Exact spin symmetry limit: S = V 

The condition for establishing the exact symmetry boundary 
implies that the solution is of the form 

(44) 

m 

So, by making S = V (or equivalently A (r) = 0 [46, 47]) in 
Eq. (10) and using the solution (44) in Eq. (43), the equation 
for fm (r) is found to be 


I{r) 


_ ^44 ^2sX+\ 

2sX +1 


Eor E = —M, we get 

Ur)=a^A'^-^\ 

gM=[b--a-H{r)\ 


(38) 


(39) 

(40) 


(f' Id m 

dP- r dr r 




\ 2 




- - + - 


-f e"^Bo(j) — es ( Bq -f (j) — ^ 


fm{r) 


+ [M^-E^+2 {E +M) y] f,n (r) = 0. 


fm{r) 

(45) 


Assuming V (r) as in Ref. [32], i.e., of the form 


where 


V(r)=ar^ + -^, 


H{r) 


2M 

—2sX +1 


—2i'.t + l 

' 7 


(41) 


Eq. (45) becomes 


Unlike from cases of Eqs. (32), (33), (34) and (35), if we 
impose that Bo and V (r) are zero, there no exist bound-state 
solutions of square-integrable. In other words, for any values 
of s, m and X in Eqs. (36)-(41), the integral (14) diverges. 


Hfm {r) = k^fm {r ), 
with 


H = EIo — estj) 


5{r) 


(46) 


(47) 


(48) 


4 Equation of motion and analysis of symmetries 


In this section, we investigate the dynamics for E ^ ±M. Eor 
this, we choose to work with Eq. (5) in its quadratic form. 
After application of the operator 


Ho 


(f- Id 2 2 


where 


(49) 


/3[(M + S(r))+/3(E-V(r)) + r-7rJ, (42) 


we get 

|p^ - 2e [(Ai -f A2) • p] -f (Ai + A2)^| \j/ (r) 

+ {[M + 5(r)]2-[£-y(r)]2-ei(j.(Bi+B2)}v/(r) 


dS (r) ,dV (r) \ / n „ 

-3^g2 + i j gi \i/{r)= 0 . 


(43) 


In this stage, it is worthwhile to mention that the Eq. (43) is 
the correct quadratic form of the Dirac equation with mini¬ 
mal, vector and scalar couplings, because the Pauli term is 
considered. 


{in — +2b{E +M), 

(50) 

^ +2a (E+M), 

(51) 

meBo — e^Bo(j)+ esBo + {E^ — ). 

(52) 


As pointed out in Ref. [32], the potential V (r) in Eq. (46) de¬ 
scribes an anharmonic oscillator. This model is a particular 
case of a proposed in Ref. [10] to study the Landau quanti¬ 
zation and the Aharonov-Bohm effect in a two-dimensional 
ring as an exactly soluble model. The model considered in 
Ref. [10] has an advantage because, besides the model con¬ 
sidered here, it also describes other physical systems, such 
as a one-dimensional ring, a straight 2D wire, a single quan¬ 
tum dot for and an isolated antidot. 
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4.2 Exact pseudospin symmetry limit: S = —V 

In this case, the condition for establishing the exact pseu¬ 
dospin symmetry limit implies that the resolution is related 
to the down component of the spinor in Eq. (11), namely 

W2 = iY.gmir)e‘^’"+^^r (53) 

m 

By making S = —V (or equivalently £ (r) =0 in Eq. (9) and 
again using Eq. (53) in Eq. (43), the equation for (r) can 
be found 


Hgm (?•) = ~k^8m (r) , 

(54) 

with 


- - 5{r) 

H = Ho es^ ^ \ 
r 

(55) 

d^ Id 9 9 

Ho = -^--^ + -2+T r^, 

dG r dr 

(56) 

where 


= {m + s — e^)^+2b{E—M), 

(57) 

172 = ^+ 2 a(£-M), 

(58) 

= [m + s) eBo — e^Bo(j) + esBo — {M^ — . 

(59) 


5 Self-adjoint extension analysis 


where Am G IR is the self-adjoint extension parameter. Eor 
Am = 0, we have the free Hamiltonian (without the 5 func¬ 
tion) with regular wave functions at the origin, and for Am ^ 
0 the boundary condition in Eq. (60) permit an singular¬ 
ity in the wave functions at the origin. 

6 The bound state energy and wave function 

In this section, we determine the energy spectrum by solving 
Eq. (47). Eor r ^0, the equation for the component /m(r) 
can be transformed by the variable change p — r]r^ resulting 
in 

P/m(p) +/m(P) - (^ + f /“(P) = 0. (63) 

Due to the boundary condition in Eq. (60), we seek for regu¬ 
lar and irregular solutions for Eq. (63). Studying the asymp¬ 
totic limits of Eq. (63) leads us to the following regular (-f) 
(irregular (—)) solution 

/m(p)=p±Je-iE(p). (64) 

With this, Eq. (63) is rewritten as 

pE"(p) + (l±v-p)E'(p)- E(p) =0. 

(65) 


In this section, we review some concepts on the self-adjoint 
extension approach. An operator with domain is 

said to be self-adjoint if and only if ^ and = 

being the adjoint of operator G. Eor smooth 
functions, ^ G C2r(lR^) with ^ (0) = 0, we should have = 
//o^, and it is possible to interpret the Hamiltonian //q (49) 
as a self-adjoint extension of //q (r 2/{0}) [48-50]. The self- 

adjoint extension approach consists, essentially, in extend¬ 
ing the domain of S!{G) in order to match ^{G"^). Erom the 
theory of symmetric operators, it is a well-known fact that 
the symmetric radial operator Hq is essentially self-adjoint 
for V > 1, while for v < 1, it admits an one-parameter family 
of self-adjoint extensions [51], where is the self- 

adjoint extension parameter. To characterize this family, we 
will use the approach in [52, 53], which is based in a bound¬ 
ary conditions at the origin. All the self-adjoint extensions 
Hq x,,, of Hq are parameterized by the boundary condition at 
the origin 




(60) 

with 



'f'o= lim r''/m(0) 

r-!-0+ 

(61) 

•El = lim ^ 
r^0+ r'' 


, (62) 


Equation (63) is of the confluent hypergeometric equation 
type 

zF”{z) + {b-z)F’{z)-aF{£)=Q. (66) 


In this manner, the general solution for Eq. (63) is 

fm{r) = F {d+,l+V,p) 

+ b,np^ie-z F{d^,l-v,p), 

with 

1 ± V 


(67) 


( 68 ) 


In Eq. (67), F{a,b,z) is the confluent hypergeometric func¬ 
tion of the first kind [45] and and h„, are, respectively, the 
coefficients of the regular and irregular solutions. 

In this point, we apply the boundary condition in Eq. 
(60). Doing this, one finds the following relation between 
the coefficients am and bm- 




Qm 



1 

■''-m''- 2-2v 

77 ^- 7 1 ™ r 

4(1-V)r^0+ 


(69) 


We note that limr^o+ diverges if v > 1. This condition 
implies that bm must be zero if v > 1 and only the regular 
solution contributes to fm{r). Eor V < 1, when the operator 
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2 3 

Bo(T) 


Equation (72) implicitly determines the bound state energy 
for the system for different values of the self-adjoint ex¬ 
tension parameter. Two limiting values for the self-adjoint 
extension parameter deserves some attention. For X,„ — 0, 
when the 5 interaction is absent, only the regular solution 
contributes for the bound state wave function. In the other 
side, for X,„ — oo, only the irregular solution contributes for 
the bound state wave function. For all other values of the 
self-adjoint extension parameter, both regular and irregular 
solutions contributes for the bound state wave function. The 
energy for the limiting values are obtained from the poles of 
gamma function, namely. 


Fig. 1 Plots of the energy {A (r) = 0) as a function of the magnetic 
field Bo for i = 1 and different values of n and m: n = 0 [solid line], 
n = 1 [dashed line] and n = 2 [dotted line]. 



Fig. 2 Plots of the energy (A (r) = 0) as a function of the magnetic 
field Bo for i = — 1 and different values of n and m: n = 0 [solid line], 
n — I [dashed line] and n = 2 [dotted line]. 


Ho is not self-adjoint, there arises a contribution of the irreg¬ 
ular solution to fm{r) [34, 54-57]. In this manner, the con¬ 
tribution of the irregular solution for system wave function 
steams from the fact that the operator Hq is not self-adjoint. 
For fm{r) be a bound state wave function, it must vanish at 
large values of r, i.e., it must be normalizable. So, from the 
asymptotic representation of the confluent hypergeometric 
function, the normalizability condition is translated in 

b„r _ r(i + v)r(r/_) 
a„, r{i-v) r{d+)' 

From Eq. (69), for v < 1, we have 

^ 3 n'' 

— // . 

dm 

Using this result into Eq. (70), one finds 

r{d+) _1 r(i-t-v) 

r(uf_) A„, 7 ^r(i-v)' 


(70) 

(71) 

(72) 


( d+ = —n for Xm = 0, (regular solution), 

[ = —n for Xm = oo, (irregular solution), 


(73) 


with n a nonnegative integer, n = 0,1,2,.... By manipulation 
of Eq. (73), we obtain 


E^-M^ = 2\ ^-^+2a{E+M) 


2n 4-1 ± {m — e(j))^ +2b {E +M) 
e^Bo^ —meBo — esBo-, S = V, 


(74) 


E^-M^ = 2\ 


/ 

- ^+2a {E-M) 


2n -f 1 ± \J{m~\-s — +2b {E —M) 


+ e^BoCl) — {m + s)eBo —esBo, S =—V. 


(75) 


As a illustration the profiles of the energy under the exact 
spin symmetry limit (S = F) as a function of the magnetic 
held Bo and with spin projection parameter s = 1 and s = — 1 
are shown in Fig. 1 and Fig. 2, respectively. From Fig. 1 and 
Fig. 2 we can note that the ground state n = 0 correspond to 
the lowest energies, as it should be for particle energy levels. 

In particular, it should be noted that for the case when 
V > 1 or when the 5 interaction is absent, only the regu¬ 
lar solution contributes for the bound state wave function 
{bm = 0), and the energy is given by Eq. (74) using the plus 
sign. The unnormalized bound state wave functions for our 
problem are 

/„,(r) = 

X F (-«, 1 ± V^, , (76) 

for S = V, and 

X F (-n, 1 ± V^, V^r^) , 


(77) 
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for S = —V, where Nf and where Ng 
stants, and 

are normalization con- 

e^B^ 

Ai = -^+2a{E+M), 

(78) 

A 2 = {m —e^)'^+2b{E+M), 

(79) 

e^B^ 

Bi = -^+2a{E-M), 

(80) 

B 2 = {m-\-s — -\-2b{E — M) 

(81) 


The self-adjoint extension is related with the presence of 
the 5 interaction. In this manner, the self-adjoint extension 
parameter must be related with the 5 interaction coupling 
constant (j)s. In fact, as shown in Refs. [34, 58] (see also 
Refs. [54, 59]), from the regularization of the 5 interaction, 
it is possible to find such a relationship. Using the regular¬ 
ization method, one obtains the following equation for the 
bound state energy 

r{d+) 1 f ^s + av \ 1 r{i + v) 

r{d-) rl^ \^s-av) Y'’r{l-v)' ^ 

By comparing Eqs. (72) and (82), this relation is found to be 


1 _ 1 / ^s + av\ 
Am rg'' \^s-avj 


(83) 


where ro is a very small radius which comes from the 5 
regularization [34, 58]. The result of Eq. (83) provides an 
explicit formula for the self-adjoint extension parameter A,,,. 
We have, therefore, derived the most important quantities for 
the system without any arbitrary parameter coming from the 
self-adjoint extension method. 


7 Nonrelativisitic limit 

Let us now examine the nonrelativistic limit of Eq. (43) by 
setting E — M + S', with S', for both cases 5 = V and 
S = —V. After applying this limit, we find 


where 


= (m — -\-AMb, 

(87) 

., 2 d 2 
•Y C tj(\ 

0 -Y4Ma, 

4 

(88) 

P = meBo — e^Boij) 4- esBo -f 2MS. 

(89) 


On the other hand, for S = —V, the term involving the poten¬ 
tial is now identically zero. The resulting equation is given 
by 


■ £/2 1 ^ 

dr^ r dr 

-<S’8m (r) =0, 




8mir) 


where 

= (m + s — e(l))^ , 

2d2 


e^B, 


0 


P = eBo {m + i) — e^Bo^ + esBo + IMS, 


(90) 

(91) 

(92) 


In order to determine the energy spectrum, we use the same 
technique above. Performing the same steps as for the rela¬ 
tivistic case, one obtains the energy levels 


= 


1 / e^Bl 


M 


^-)-4Ma 


2n -f 1 ± \J{m — e<^ f'+AMb 


+ ^ -{m + s) eBo - esBo] , S = V, 


(93) 


S = (2n -I-1 ± \m - 1) 

+ ^ [e^BoCj) - {m + s) eBo - esBo] , S =-V, 
The corresponding wave functions are given by 


fm{r) = 


Pb2 


-|-4Ma 


X r 


.ip {fn—e(l))^+4Mb 


±|-v/ {m—e^)^+4Mb 


^+4Mar^ 


(94) 


Hxff = IMSxff, 
where 

H = {p — eA)^ — esa ■ B + 2M [5 (r) + V (r)]. 


(84) 


(85) 



for S = V, 


Using the ansatz of Eq. (11) in Eq. (84), again, we get the 
equation for /,„ (r) (for S = V) 


d^ 

dr^ 


Id -2 2 A ^('■) 

- — -f +1] ^ 

r dr r 


fm{r) 


-^/m(r)=0. 


( 86 ) 


fm{r) = 


r2 
e Bn. 


n ± I |m— 


+ \m-eij)\ -^eBor^ 


X F ( —n, 1 ± |m — 1, -eBor ) , 


for S = -V. 


(96) 
































8 Conclusions 6. 

In this paper, we have studied the relativistic quantum dy- 7. 
namics of a spin-1/2 charged particle with minimal, vector 
and scalar couplings. The minimal coupling was chosen as 8. 
being one which leads to the spin-1/2 Aharonov-Bohm ef¬ 
fect. In a first attempt, we have solved the equation of first 9. 
order Dirac. We verified that there are isolated solutions for 
the system for some special cases. These solutions depend 10. 
on the values assumed by the spin projection parameter s, 
as well as on the choice of the scalar S (r) and vector V (r) 11. 

potential functions. 

In contrast to what was addressed in the literature, we 12. 
have considered the correct quadratic form of the Dirac equa¬ 
tion with minimal, vector and scalar couplings. As we have 13. 
mentioned before, in approach of Ref. [32], the authors have 
not taken into account the term that depends explicitly on 
the spin in the Pauli equation of motion. We have revisited 14. 
the dynamics of the system in detail and show that the cor¬ 
rect approach should involve the spin element. In this order, 15. 
we have derived the Pauli equation and studied the dynam¬ 
ics of the system taking into account exact symmetry spin 
and pseudospin limits. Because of the equation of motion in- 16. 
eludes a d function, we have used the self-adjoint extension 
method to specify the proper boundary conditions at the ori- 17. 
gin. The analytical solutions of the model allow a calculation 
of expressions for the energy eigenvalues and wave func¬ 
tions for both symmetry limits. We verify that the presence 18. 
of the spin element in the equation of motion introduces a 
correction in the expressions for the bound state energy and 
wave functions, a fact that does not occur in Ref. [32]. 19. 
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